Entanglement and Bell States in Superconducting Flux Qubits 
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We theoretically study macroscopic quantum entanglement in two superconducting flux qubits. 
To manipulate the state of two flux qubits, a Josephson junction is introduced in the connecting loop 
coupling the qubits. Increasing the coupling energy of the Josephson junction makes it possible to 
achieve relatively strong coupling between the qubits, causing two-qubit tunneling processes to be 
even dominant over the single-qubit tunneling processes in the states of two qubits. It is shown that 
due to the two-qubit tunneling processes both the ground state and excited states of the coupled 
flux qubits can be a Bell type state, maximally entangled, in experimentally accessible regimes. The 
parameter regimes for the Bell states are discussed in terms of magnetic flux and Josephson coupling 
energies. 
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PACS numbers: 74.50.+r, 85.25.Cp, 03.67.-a 

Introduction.— Quantum entanglement is a fundamen- 
tal resource for quantum information processing and 
quantum computing Numerous proposals have been 
made for the creation of entanglement in solid-state sys- 
tems: quantum dots 0, 0, 4, 5], Kondo impurities p], 
carbon nanotubes [3] and so on. For superconducting 
qubits coherent manipulation of quantum states in a con- 
trollable manner has enabled to generate partially entan- 
gled states @, 1^, [13] ■ However, in realizing such quan- 
tum technologies, highly entangled quantum states such 
as the Bell states of two qubits are required Of par- 
ticular importance, therefore, are the preparation and 
measurement of such maximally entangled states. 

Recently, several types of superconducting qubits have 
been demonstrated experimentally. In superconducting 
qubits experimental generations of entanglements have 
been reported for coupled charge qubits |8j] and coupled 
phase qubits , but the maximally entangled Bell states 
is far from experimental realization as yet. This pa- 
per aims to answer on how maximally entangled states 
can be prepared by manipulating system parameters in 
solid-state qubits. We quantify quantum entanglement 
in the two superconducting flux qubits based on the 
phase-coupling scheme. In this study we show that, if 
the coupling strength between two flux qubits is strong 
enough, simultaneous two-qubit coherent tunneling pro- 
cesses make it possible to create a maximally entangled 
state in the ground and excited states. Furthermore, it is 
shown that the ranges of system parameters for a max- 
imally entangled state are sufficiently wide that a Bell- 
type state should be realizable experimentally. Actually, 
a coherent two-qubit flipping processes has been experi- 
mentally observed in inductively coupled flux qubits [13] • 
However, the strength of the inductive coupling ll| is too 
weak to achieve a maximally entangled state. 

We use a phase-coupling scheme [l3] to obtain suf- 
ficiently strong coupling and show a maximally entan- 
gled state between two flux qubits. Very recently, this 



phase-coupling scheme has been realized in an experi- 
ment [isj using two four-j unctions flux qubits. To the- 
oretically study a controllable coupling manner in the 
flux qubits, the phase coupling scheme has also been 
employed [ij, [l5|. Further, there have been stud- 
ies about somewhat different types of phase-coupling 
schemes [13, [13 ■ The phase-coupling scheme for two flux 
qubits (See Fig. [T]) is to introduce a connecting loop in- 
terrupted by a Josephson junction in order to couple two 
three-junctions qubits [3, [l3]- In the connecting loop, 
the Josephson energy depends on the phase difference, 
(f'Li — 4>Ri^ and the coupling energy, E'j, of the junc- 
tion. As varies E'j, the coupling strength between the 
flux qubits, defined by the energy difference between the 
same direction current state and the different direction 
current state in the two flux qubits, can be increased to 
be relatively strong. The coupling strength between the 



X2 



Ej 
Ej 



1 rt§n_ 



'm 



X3 



FIG. 1: Left and right superconducting loops are connected 
each other by a connecting loop interrupted by a Josephson 
junction. The state of each qubit loop is the superposed state 
of the diamagnetic and paramagnetic current states assigned 
by II) and |t), respectively, which can make the loop being 
regarded as a qubit. For example, the state, jit), out of four 
possible basis of two-qubit current states is shown, where the 
arrows indicate the flow of Cooper pairs and thus in reverse 
direction is the current. Here, (oppositely ®) denote the 
directions of the magnetic fields, = ^l{r)/^o, in the 

qubit loops. Ej\, Ej, and E'j are the Josephson coupling 
energies of the Josephson junctions in the qubit loops and 
the connecting loop and iji's are phase differences across the 
Josephson junctions. 
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phase-coupled flux qubits is a monotonously increasing 
function of E'j [12,]. 

For small Ej{<^ Ej), the effective potential for the 
two qubits has a symmetric form in the phase variable 
space. Then single-qubit tunneling processes in the two 
qubit states are predominant in the effective potential. 
They have almost the same tunneling amplitudes. The 
eigenstates are nearly degenerate at the operating point 
of the external fluxes so that the entanglement between 
two qubits is very weak. As increases the Josephson cou- 
pling energy, Ej, the shape of the effective potential is de- 
formed to be a less symmetric form. Then the amplitudes 
of single-qubit tunneling processes are decreased. How- 
ever, the deformation of the effective potential allows a 
two-qubit tunneling process not negligible and even dom- 
inant over the single-qubit tunneling processes. We will 
then derive explicitly the contribution of two-qubit tun- 
neling processes. To generate a Bell type of maximally 
entangled state, the two-qubit tunneling processes are 
shown to play an important role. 

Model.— Let us start with the charging energy of the 
Josephson junctions: Ec = ^ ($o/27r)^ (Z]Li(^ii'^L + 
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CBi4>Ri) + C"0' ), where CL(R)i and C" are the capaci- 
tance of the Josephson junctions for the left (right) qubit 
loop and the connecting loop, respectively. Here, (/)'s are 
the phase differences across the Josephson junctions and 
$0 = h/2e the superconducting unit flux quantum. Since 
the number of excess Cooper pair charges on the Joseph- 
son junctions, Ni = Qi/2e with Qi = C(<i>o/27r)(/)i, is 
conjugate to the phase difference (pi such as [^i, Ni] — i, 
the canonical momentum, Pi = Nih = —ihd/d4>i, can be 
introduced. Then the Hamiltonian is given by 

H = \PlM-^'p, + UM4>), (1) 

where = (<i>o/27r)^Ci5ij is the effective mass and 

4> = {(t>Li,4'L2,(l>L3,(l)Ri,4'R2j<l)m, (!>')■ If we ueglect the 
small inductive energy, the effective potential becomes 
the energy of the Josephson junctions such as Uce{4>) — 
Eti EMI - cos^lO +ELi EMI - cos^ffi) + E'j{l - 
cos 0'). 

Since the qubit operations are performed experimen- 
tally at near the co-resonance point, we can set = 
fji = 0.5 with = $L(fl)/$o and the total flux 

^L{R) threading the left (right) qubit loop. In exper- 
iments, the two Josephson junctions with phase differ- 
ences, 4'L{R)2 and (j)L(R)3, are considered nominally the 
same so that it is reasonable to set Ej2 = Ej3 = Ej and 
4>L{R)2 — 4>L{R)z- Furthermore, if one can neglect the 
small inductive flux, the boundary conditions in the left 
and right qubit loops and the connecting loop are given 
approximately, 2iT{nL(R) + fL(R)) - {(I)l(R)i + 4>l(r)2 + 
4'L(R)7i) — and 27rr + {4>li — (t'Ri) —(!)' = with inte- 
gers uljUr and r. Introducing the rotated coordinates 
with (j>p = (0L3 + 0fl3)/2 and 0„ = (0^3 - cj)R,z)/2, 




FIG. 2: (Color online) The effective potential of Eq. ^ 
for the coupled flux qubits at the co-resonance point /l = 
/h = 0.5 as a function of <^l3 and (jiRi with Eji = Q.lEj. 
In the contour plots of the effective potential, the four local 
minima correspond to the current states of the coupled flux 
qubits: ITT) , ITi) , lit), and The left-right arrows (^) 
on the effective potential proflles indicate some of single-/two- 
qubit tunneling processes (For example, ||T) ^^=^ ITT) and 
ITT) ITT)) respectively), (a) For E'j — 0, the four local 

minima of the effective potential have nearly the same ener- 
gies. The potential barrier of the double well for a two-qubit 
tunnehng process are much wider and higher than the poten- 
tial barriers for the single-qubit tunneling processes. Thus, 
single-qubit tunneling processes are dominant and the entan- 
glement is very weak, (b) For E'j = Q.GEj, only |TT) and 
ITT) are lifted up to a higher energy while |TT) and |TT) stay 
at the same energies. Then the potential barrier for the two- 
qubit tunneling process between |Ti) and |TT) is not changed 
so much that the tunneling amplitude also is not changed. 
An asymmetric double well potential for single-qubit tunnel- 
ing processes makes the tunneling amplitude decreasing as E'j 
increases. Therefore, the two-qubit tunneling process plays a 
significant role on the entanglement of the coupled fiux qubits. 

the Hamiltonian can be rewritten in the form, H — 
Pp/2Mp + P^/2M„ + Ucgi(t>p,4>m): where the effective 
potential becomes 

C/cff(0p,0m) ^ 2Eji{l + COS2(f>pCOS2(f>jn) (2) 

+AEj{l — COS0P cos (/)„,.) + Ej{l — cos4(/)„j). 

Here, Mp = h^{2Ci+C)/e^, M„ = h^{2Ci + C + 4C')/e^ 
and the conjugate momentum is Pp(m) = —ihd / d(/)p(^jn) 
with commutation relation, [4>p{m) , Pp(,n)] = i^- 

We display the effective potential as a function of 
and (j)R3 in Fig. [21 The effective potential is shown to 
have the four local potential minima corresponding to 
the four states of the coupled qubits, i.e., |tT) > ITi) j lit)) 
and III). At the local minima of the effective potential, 
one can obtain the energy levels: Ess' = {^/'2)fi{'-^p,ss' + 
^rn.ss') + Ueei4'p,ss',4>m,ss'), whcrc s,s' = 1/2 and -1/2 
stand for spin up, |t), and spin down, ||), respec- 
tively. In the harmonic oscillator approximation (20| , the 
characteristic oscillating frequencies, (jJp(m),ss', are given 

by UJp(^rn),ss' = {kp(m),ss' IMp{^rn)Y/'^ with kp(m),ss' = 
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52C/efi(0p,(^m)/5(^p(„)|(0^_^^,^0^^^,), where 0p,ss'(0m,ss') 
are the values of (t'p{4'm) at the local potential minimum. 
Then the tight-binding approximation gives the Hamil- 
tonian in the basis, {|TT) , |Ti) . lit) , lii)}, as follows; 



H=J2 E,s'\s,s'){s,s'\+Hl + H^, (3) 

s,s'—±^ 

Hl=Y. {- ■'^'1 - *f 1^' -^'l) ' (4) 

^2= E {-t^2l'^,''^){^S-s\-4\s,-s){-S,s\),i5) 

where H^^^ describe the single(two)-qubit tunneling pro- 
cesses between the two-qubit states with the tunneling 
amplitudes ^{''^(ij''')- Normally, the single-qubit tunnel- 
ing amplitudes are much larger than the two-qubit tun- 
neling amplitudes, i.e., ij'** t^'^ ■ Thus, the two-qubit 
tunneling Hamiltonian _ff J can be neglected when Ej=0. 

Two-qubit tunnelings and entanglements.— For the 
weak coupling limit, Ej ^ Ej, the single-qubit tunnel- 
ings between the two-qubit states can be obtained by the 
tunnelings in the well-behaved double well potentials as 
shown as left-right arrows in Fig. [2ta). As increases Ej, 
the well-behaved double well potentials for the single- 
qubit tunnelings become an asymmetric double well po- 
tentials shown clearly in Fig. [H^b). Then the asymmetry 
of the double well potentials makes the single-qubit tun- 
neling amplitudes decrease drastically. However, as the 
blue line shown in Fig. [2{b), the corresponding dou- 
ble well potential to the two-qubit tunneling between the 
states III) and |||) is not changed much qualitatively. 
Actually, from the effective potential in Eq. ^ , the dou- 
ble well potential is given by 

C/cff(</'p,0) = 2i;ji(l+cos20p)+4i;j(l-cos(/)p). (6) 

The two-qubit tunneling amplitudes if between the 
states Itt) and |||) is written in the WKB approximation 
f20l. 



ZTT 



exp 



E, 



(7) 



Here Mp, ujp^ss and Ucff{(f>p, 0) do not depend on the ca- 
pacitance, C", and the Josephson coupling energy, E'j, 
of the Josephson junction in the connecting loop. Equa- 
tion ([7|) shows that the two-qubit tunneling amplitude 
remains unchanged even though Ej is increased. As a 
consequence, the two-qubit tunneling between the states 
III) and III) in Hamiltonian _ffj of Eq. ([5]) will play 
a crucial role to improve the entanglement between the 
coupled flux qubits. 

For Ej — 0, we obtain the tunneling amplitudes; 
tf '-^ « 0.0075£;,7 and = « 0.00024£;,7 with Eji = 
O.lEj. As expected, t'^''' < t^'^^ and the two-qubit tun- 
neling terms of the Hamiltonian are negligible. When 



E'j ^ 0, in order to get the single-qubit tunneling ampli- 
tudes between an asymmetric double well potentials, the 
Fourier Grid Hamiltonian Method 2l| is employed be- 
cause the WKB approximation cannot be applicable in 
the asymmetric double well potential. For Ej = 0.6Ej 
as shown in Fig. [Sfb), we obtain the tunneling ampli- 
IQ-^Ej and « 0.00024£;j. Here, the 



tudes of ti'^ 



tunneling amplitude is neglected because the the wave 
function overlap between the states |t|) and ||t) is neg- 
ligible in Fig. lib). 

One measure of entanglement is the concurrence, C, 
for an arbitrary state of coupled two qubits. The con- 
currence ranges from for nonentangled to 1 for max- 
imally entangled states. For a normalized pure state, 
IV") = a|ii) + &|Ti) + c||T) +of|TT), the concurrence [H 
is given by 



C{\ip))^2\ad-bc\. 



(8) 



We evaluate the concurrence in Eq. ([5]) numerically to 
show that a maximally entangled state is possible in the 
ground state by varying the Josephson coupling energy 
in the connecting loop. In Fig. [2Ia), we plot the concur- 
rence for the ground state, C{\iPg)), of the Hamiltonian 
in Eqs. ([S])-® as a function of fpi^m) = (/l ± /i?)/2 
for E'j = 0.6Ej and Eji = 0.7Ej. A broad ridge of 
the concurrence along the line fp = 0.5 shows a high 
entanglement between the two qubits. This region corre- 
sponds to the central part of the honeycomb type poten- 
tial of coupled qubits in Ref. [l^ near the co-resonance 
point, /l = fa = 0.5. Away from the co-resonance point 
we can numerically calculate concurrences without using 
the analytic effective potential in Eq. ^ obtained by in- 
troducing a few approximations [l2| . Figure [3Jb) shows 
the cut view of the concurrence at fp = 0.5 (dotted line 
in Fig. El^a)) as a function of /,„ for various values of 
E'j. For E'j = 0.005£'j, corresponding to the couphng 
strength of inductively coupled qubits in the experiment 
of Ref. X^i a partial entanglement can only exist around 
fm ~ 0. As increases E'j, i.e., when the two-qubit tun- 
neling becomes dominant over the single-qubit tunneling, 
a maximum entanglement appears. In Fig. ^c), we dis- 
play the coefficients of the ground state wavefunction, 

IV'g) = aifm) III) + bifra) |T|) + c(/„) ||T) + d(/„0 |TT), 

as a function of for E'j = 0.2Ej. We see the ground 
state wavefunction: 



for U < -0.04 




(Ui) + ITT)) for -0.04 </„< 0.04 . (9) 
for /,„ > 0.04 



This ground state has the Bell type of maximally entan- 
gled state, |$+) = (III) + |TT))/\/2, for -0.04 < fm < 
0.04. Similarly, one can see that the 1st excited state 
shows, actually, another Bell type of maximally entan- 
gled state, 1$') = (|||) - |TT))/\/2. 
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FIG. 3; (a) Concurrence for the ground state of the coupled 
flux qubits as a function of fp(m) = (/i ± fR)/'^ for E'j — 
O.GEj and Eji = Q.7Ej. A high entanglement between two 
flux qubits is seen in the broad ridge with /p = 0.5. In the 
contour plot, ITi), UT), and denote the coupled- 

qubit states located far from the co-resonance point, /l — 
/ij = 0.5. (b) Concurrences for various E'j along the line of 
fp = 0.5 (dotted line in (a)) for Eji = 0.7Ej. As increases 
e'j, the maximum entanglement appears for E'j > 0.03 -Ej. 
(c) The coefficients of the ground state wavefunction, lipc) ~ 
aUm) lU) +&(/-) ITi) +c(/™) UT> |TT>, as a function 

of /„. for the case E'j = 0.2Ej in (b). For -0.04 < < 0.04, 
the Bell type of maximally entangled state, \^^) ~ + 
|tt))/\/2, is shown clearly because a = d = l/\/2 and h = 
c = 0. (d) Concurrences for various Eji along the line of 
fm = Q for E'j = 0.6-Ej (dashed line in (a)). As decreases 
Ej\, since the two-qubit tunneling amplitude increases, the 
width of concurrence peak becomes so wide that the maximal 
entanglement should be observable experimentally. 



coupling between the qubits and if''^ ^ t2'^. As E'j in- 
creases, the coupling between the qubits is much stronger 
and ^2 becomes much dominant over than other tun- 
neling processes, if, if, and t\. In a strong couphng 
limit, tf,tf,t2 < ^2 and Eg^s < Es-s- If we with 
the split energies, E^ — Eq — 6E and = Eo + 6E 
around the co- resonance point, /l — fn ~ 0.5 Then the 
concurrence of the ground state is given by C{\^l>G)) — 
\t%l \J (^2)^ + {^EY\^ which corresponds to the behaviors 
in Fig. [^d). At the co-resonance point (/l = Jr = 0.5), 
the ground state is in a maximally entangled state, i.e., 
C(|'0g)) = 1- This shows the important role of the two- 
qubit tunneling process making the two qubits entangled 
and maintaining the highly entangled state against fluc- 
tuations of Jl^h away from the co-resonance point as long 
as SE < q. 

Summary.— We studied the entanglement to achieve a 
maximally entangled state in a coupled superconducting 
flux qubits. A Josephson junction in the connecting loop 
coupling the two qubits was employed to manipulate the 
qubit states. As increases the Josephson coupling en- 
ergy of the Josephson junction, the two-qubit tunneling 
processes between the current states play an important 
role to make the two flux qubit strongly entangled. It was 
shown that a Bell type of maximally entangled states can 
be realized in the ground and excited state of the coupled 
qubit system. We also identified the system parameter 
regime for the maximally entangled states. 
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The cut view of the concurrence at = (dashed 
line in (a)) is shown in Fig. ^d) for various Eji with 
E'j = 0.6Ej. As Eji increases, the barrier of the dou- 
ble well potential in Eq. ([6]) becomes higher. Thus the 
two-qubit tunneling amplitude, decreases such that 
q/Ej w 0.0016, 0.0008, 0.00024 for Eji/Ej^O.65, 0.67, 
0.7, respectively. In Fig. E^d), as a result, the width of 
the concurrence peak becomes narrower, which will make 
an experimental implementation of the maximum entan- 
glement difficult. On the other hand, if the value of Eji 
becomes much smaller than Eji — 0.65Ej, the barrier of 
double well potential can be lower than the ground energy 
of the harmonic potential wells, E^s, and the two-qubit 
states, III) and |tT), in Fig- [2Ib) will not be stable. 
In addition, as seen in Fig. ^h), the range of E'j for 
maximal entanglement is approximately E'j > 0.03Ej. 
Therefore, for a maximally entangled state of the two 
flux qubits, it is required that Eji should be controlled 
around 0.7£;j with E'j > 0.03 Ej. 

For E'j <C Ej, the concurrence is very small because 
the four states are nearly degenerated due to the weak 
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